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Abstract. Malliavin Calculus is about Sobolev-type regularity of functionals 
on Wiener space, the main example being the Ito map obtained by solving sto- 
chastic differential equations. Rough path analysis is about strong regularity of 
solution to (possibly stochastic) differential equations. We combine arguments 
of both theories and discuss existence of a density for solutions to stochastic 
differential equations driven by a general class of non-degenerate Gaussian pro- 
cesses, including processes with sample path regularity worse than Brownian 
motion. 



1. Introduction 

It is basic question in probability theory whether a given stochastic process 
{Yj : t > 0} with values in some Euclidean space admits, at fixed positive times, 
a density with respect to Lebesgue measure. In a non-degenerate Markovian set- 
ting - cllipticity of the generator - an affirmative answer can be given using Weyl's 
lemma as discussed in McKean's classical 1969 text [53]. Around the same time, 
Hormander's seminal work on hypoelliptic partial differential operators enabled 
probabilists to obtain criteria for existence (and smoothness) for densities of cer- 
tain degenerate diffusions. This dependence on the theory of partial differential 
equations was removed when P. Malliavin devised a purely probabilistic approach, 
perfectly adapted to prove existence and smoothness of densities. 

We recall some key ingredients of Malliavin's machinery, known as Malliavin 
Calculus or stochastic calculus of variations. Most of it can formulated in the setting 
of an abstract Wiener spaces (W,Tl, /i). The concept is standard [3j [22J, [HI HE] as 
is the notion of a weakly non- degenerate Revalued Wiener functional <p which has 
the desirable property that the image measure ip^fi is absolutely continuous with 
respect to Lebesgue measure on IR e . (Functionals which are non-degenerate have a 
smooth density.) Precise definitions are given later on in the text. 

Given these abstract tools, we turn to the standard Wiener space Co ([0, 1] , R d ) 
equipped with Wiener measure. From Ito's theory, we can realize diffusion processes 
by solving the stochastic differential equation 

d 

dY = Vi (Y) o dB l + V (Y) dt = V (Y) odB + V (Y) dt, Y (0) = y e K e 

i=l 

driven by a d-dimensional Brownian motion B along sufficiently well-behaved (drift 
- resp. diffusion) vector fields Vq, Vd- The Ito-map B i— » Y is notorious for its 
lack of strong regularity properties which rules out the use of any Frechet-calculus. 
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On the positive side, it is smooth in a Sobolev sense on Wiener space (" smooth in 
the sense of Malliavin"). Under condition 

(E) : span^,...,^ = T yo K e ^ K e ; 

or a less restrictive Hormander's condition (allowing all Lie brackets of Vq, ...,Vd 
in spanning 7^ R e ) the solution map B i— > Y tl for t > 0, is non-degenerate and 
one reaches the desired conclusion that Y t has a (smooth) density. This line of 
reasoning due to P. Malliavin provides a direct probabilistic approach to the study 
of transition densities and has found applications from stochastic fluid dynamics 
to interest rate theory. It also shows that the Markovian structure is not essential 
and one can, for instance, adapt these ideas to study densities of Ito-diffusions as 
was done by Kusuoka and Stroock in [T?] . 

Our interest lies in stochastic differential equations of type 

(1.1) dY = V (Y) dX + V Q (Y) dt 

where X is a multi-dimensional Gaussian process. Such differential equations 
arise, for instance, in financial mathematics [H [TOl [TT] or as model for studying 
crgodic properties of non-Markovian systems |12j . Assuming momentarily enough 
sample path regularity so that makes sense path-by-path by Riemann-Stieltjes 
(or Young) integration the question whether or not Yt admits a density is important 
and far from obvious. To the best of our knowledge, all results in that direction are 
restricted to fractional Brownian motion with Hurst parameter H > 1/2. Existence 
of a density for Y t , t > 0, was established in [55] under condition (E). Smoothness 
results then appeared in [26], and under Hormander's condition in [1]. 

The purpose of this paper is to give a first demonstration of the powerful interplay 
between Malliavin calculus and rough path analysis. After some remarks on Ti.- 
differentiability and a representation of the Malliavin covariance in terms of a 2D 
Young integral, we show that weak non-degeneracy (in the sense of Malliavin) of 
solutions to p.l[) at times T > 0, and hence existence of a density, remains valid in 
an almost generic sense. Our assumptions are 

• The vector fields V = (Vi, ■■■,Vd) at yo span T Vo M., i.e. condition (E). 

• The continuous, centered Gaussian driving signal X is such that the sto- 
chastic differential equation (jl.ip makes sense as rough differential equation 
(RDE), [18l[2T]. Applied in our context, this represents a unified framework 
which covers at once Gaussian signals with nice sample paths (such as fBM 
with H > 1/2), Brownian motion, Gaussian (semi-) martingales, and last 
not least Gaussian signals with sample path regularity worse than Brownian 
motion provided there exists a sufficiently regular stochastic arecQ- 

• The Gaussian driving signal is sufficiently non-degenerate which is clearly 
needed to rule our examples such as X = or the Brownian bridge X t — 
B t (u) - (t/T) Br (u). 

Smoothness of densities remains an open problem, some technical remarks about 
the difficulties involved are found in the last section. 

Acknowledgement 1. We would like to thank James Norris for related discus- 
sions and an anonymous referee whose comments prompted a clarification in our 
discussion of TL-regularity. The second author is partially supported by a Lever- 
hulme Research Fellowship. 



For orientation, fractional Brownian motion is covered for any H > 1/3. 
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2. Preliminaries on ODE and RDEs 

2.1. Controlled ordinary differential equations. Consider the ordinary differ- 
ential equations, driven by a smooth Revalued signal f = f (t) along sufficiently 
smooth and bounded vector fields V — (Vi, Vd) and a drift vector field Vb 

dy = V(y)df + V (y) dt, y (t„) - y Q e R e . 

We call U{^ t (yo) = yt the associated flow. Let J denote the Jacobian of U. It 
satisfies the ODE obtain by formal differentiation w.r.t. yo. More specifically, 



^- £ ULt a (yo + ea) 



£ = 



is a linear map from R e R e and we let j/<_ to (yo) denote the corresponding e x e 



-to 

matrix. It is immediate to sec that 



pU (yo) 



dt 

where ■ denotes matrix multiplication and 



d Mi (uL tQ (y ),t 



■ JLta (yo) 



± M f(y,t)=Y l vay)j t fi + Vo'(y). 

Note that J^<_ t = j( 2 ^t 1 ' Jti*-t - We can consider Gateaux derivatives in the 
driving signal and define 

Dull f -f — U f+£h 

One sees that DhU/^ satisfies a linear ODE and the variation of constants formula 
leads to 



D h u t t (yo) = f E J Ls (Vi (uLo 



dhl. 



2.2. Rough differential equations. Let p 6 (2,3). A geometric p-rough path x 
over R d is a continuous path on [0, T] with values in G 2 (R d ) , the step-2 nilpotent 
group over R d , of finite p- variation relative to d, the Carnot-Caratheodory metric 
on G 2 (R d ), i.e. 

sup sup V]d (x ti ,x ti+1 ) p < co. 

n ti<...<t„ 

Following 0101 we realize G 2 (R d ) as the set of all (a, 6) 6 R d © R dxd for 
which Sym (b) = a® 2 1 2. (This point of view is natural: a smooth Revalued path 
x = (x\) ,, enhanced with its iterated integrals J Q Jq dx^dxl, gives canonically 

rise to a G 2 (R d )-valued path.) Given (a, 6) G G 2 (R d ) one gets rid of the redundant 
Sym (b) by (a, 6) h-» (a, 6 — a®/2) S R d © so (d). Applied to a; enhanced with its 
iterated integrals over [0, t] this amounts to look at the path x and its (signed) 

areas / ( 

x l ~ x h) dx J s — Jo (%i ~ x ~oj dXg, i,j £ {1, . . . , d} . Without going into too 

much detail, the group structure on G 2 (R d ) can be identified with the (truncated) 
tensor multiplication and is relevant as it allows to relate algebraically the path 
and area increments over adjacent intervals; the mapping (a, b) <— * (a, b — a®/2) 
maps the Lie group G 2 (R d ) to its Lie algebra; at last, the Carnot-Caratheodory 
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metric is defined intrinsically as (left-)invariant metric on G 2 (M d ) and satisfies 

H + l&| 1/2 <rf((0,0),(a,6))<|a| + |6| 1/2 . 

One can then think of a geometric p- rough path x as a path x : [0, T] — > W l en- 
hanced with its iterated integrals (equivalently: area integrals) although the later 
need not make classical sense. For instance, almost every joint realization of Brow- 
nian motion and Levy's area process is a geometric p- rough path [3TJ O [5] . The 
Lyons theory of rough paths [201 HH E] then gives deterministic meaning to the 
rough differential equation (RDE) 

dy = V (y) dx 

for Lip 7 -vector fields (in the sense of Stein), 7 > p. By considering the space-time 
rough path x = (t, x) and V = (Vo, V%, Vd) this form is general enough to cover 
differential equations with drift0. The solution induces a flow yo <— * U^ to (yo)- The 
Jacobian J t ^_ to of the flow exists and satisfies a linear RDE, as does the directional 
derivative 

for a smooth path h. If x arises from a smooth path x together with its iterated 
integrals the translated rough path T/ t x is nothing but x + h together with its 
iterated integrals. In the general case, we assume h s (79-var w jth l/f>+ 1/<Z > 
1 , the translation T^x can be written in terms of x and cross- integrals between 
7Ti (xo,.) =: x and the perturbation h. (These integrals are well-defined Young- 
integrals.) 

Proposition 1. LefX be a geometric p-rough paths overR d and h £ C q ([0, 1] ,M. d ) 
such that 1/p + 1/q > 1. Then 

D h UH (yo) = f E J * * 0)) dh* 

where the right hand side is well-defined as Young integral. 

Proof. At least when 7 > p+ 1, both J^_ and D^U^q satisfy (jointly with U^_ ) 
a RDE driven by X. This is an application of Lyons' limit theorem and discussed 
in detail in [THl [20] • A little care is needed since the resulting vector fields now have 
linear growth. It suffices to rule out explosion so that the problem can be localized. 
The needed remark is that J^_ also satisfy a linear RDE of form 

rfj t l = dM x ■ J* (y ) 

where dM x = V (U£_ (y )) dX. Explosion can be ruled out by direct iterative 
expansion and estimates of the Einstein sum as in [18] . □ 

3. RDES DRIVEN BY GAUSSIAN SIGNALS 

We consider a continuous, centered Gaussian process with independent compo- 
nents X = (X 1 , X d ) started at zero. This gives rise to an abstract Wiener 
space (W,H,/J,) where W = H c C ([0,T] ,R d ). Note that H = ®f =1 H® and 
recall that elements of H are of form ht = K(Xt£(h)) where £ (h) is a Gaussian 

If x lifts a Revalued path x, then x is construct with cross-integrals of type f xdt, f tdx 
all of which are canonically denned. (2) Including Vo in the collection V leads to suboptimal 
regularity requirements for Vo which could be avoided by a direct analysis. 
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random variable. The ("reproducing kernel") Hilbcrt-structure on TL is given by 
(h,h') H :=E(t(h)t(h')). 

Existence of a Gaussian geometric p-rough path above X is tantamount to the 
existence of certain Levy area integrals. From the point of view of Stieltjes integra- 
tion the existence of Levy's area for Brownian motion is a miracle. In fact, a subtle 
cancellation due to orthogonality of increments of Brownian motion is responsible 
for convergence and this suggests that processes with sufficiently fast decorrela- 
tion of their increments will also give rise to a stochastic Levy area. The resulting 
technical conditions appears in [5T] for instance. For Gaussian processes, a cleaner 
(and slightly weaker) condition can be given in terms of the 2D variation prop- 
erties of the covariance function R(s,t) = E(X s ®X t ) = diag The 
assumption, writing X t ,v := X t > — X t , 

\ R \Um[o,^--= su p Ei E (^i.*^.^)r<~ 

D=(U) ij 

for p < 2 is known [8] to be sufficient (and essentially necessary) for the existence 
of a natural lift of X to a geometric p-rough path X for any p > 2p. Observe that 
the covariance of Brownian motion has finite p- variation with p = 1. As a more 
general example, a direct computation shows that fractional Brownian motion has 
finite p-variation with p = 1/ (2H). 

The assumption of l-Rlp.va,. < oo has other benefits, notably the following em- 
bedding theorem [8]. Since it is crucial for our purposes we repeat the short proof, 
as we apply it componentwise we can assume d = 1. 

Proposition 2. Let R be the covariance of a real-valued centered Gaussian process. 
If R is of finite p-variation, then TL °-> C p ~ var . More, precisely, for all h £ TL, 



\ h \p-var:[s,t] ^ \J ( h , h )n\J ' R p-var;[s,t] 2 



Proof. Every element h £ TL can be written as ht — E (ZXt) for the Gaussian r.v. 
Z = £ (h). We may assume that (h, h) n = E \ Zr) = 1. Let (tj) be a subdivision of 

[s, t] and write \x\, r = x i) f° r r > 1. If p' denote the Holder conjugate of p 
we have 

E l^i+i r = su p E^-^.^+i = sup E rE^A^+i 



< sup fe2l3jl3 k E(X tj , tj+1 X t , 



™p \ lEi^-r'^n Ei E (^^ +1 ^^ +1 )i 

/ \ 1/( 2 P) 

^ ( E \ E ( x t 3 ,t J+1 X tk , tk+1 )\ P J < y / i? p _ var;[Sit] 2. 
Optimizing over all subdivision (tj) of [s,t] , we obtain our result. □ 
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One observes that for Brownian motion (p = 1) this embedding is sharp. Fur- 
thermore, the rough path translation T^X, which involves Young integrals, makes 
sense if p < p* —3/2. (This critical value comes from 1/p+l/ (2p) ~ 1/p+l/ (2p) = 
3/ (2p) and equating to 1.) The moral is that one can take deterministic directional 
derivatives in Cameron-Martin directions as long as p < 3/2. (En passant, we 
see that the effective tangent space to Gaussian RDE solutions is strictly bigger 
than the usual Cameron-Martin space as long as p < 3/2. In the special case of 
Stratonovich SDEs a related result predates rough path theory and goes back to 
Kusuoka [IS]). For this reason we will assume 

p < 3/2 

for the remainder of this paper. This entails that we are dealing with geometric 
p-rough paths X for which we may assume 

pe (2/0,3). 

Definition 1. [13], [531 Section 4.1.3], [30] Section 3.3]Given an abstract Wiener 
space (W,TL,p,), a random variable (i.e. measurable map) F : W — > K is continu- 
ously TL- differentiable, in symbols F G C^, if for p-almost every u>, the map 

h G TL >-> F [u + h) 

is continuously Frechet differentiable. A vector-valued r.v. F — [F 1 , ...,F e ) : W — > 
IR e is continuously TL- differentiable iff each F l is continuously TL-differentiable. In 
particular, p-almost surely, DF (iS) — (DF 1 (w) , ...,DF e (a;) J is a linear bounded 
map from TL — > R e 

Remark 1. (1) The notion of continuous TL- differentiability was introduced in [13] 
and plays a fundamental role in the study of transformation of measure on Wiener 
space. Integrability properties of F and DF aside, C^-regularity is stronger than 
Malliavin differentiability in the usual sense. Indeed, by |24[ Thm 4.1.3] (see also 
[13j . [301 Section 3.3],) implies D ^-regularity where the definition of ID^ C is 
based on the commonly used Shigekawa Sobolev space D 1,p . (Our notation here fol- 
lows [241 Sec. 1.2, 1.3.4]/ This remark will be important to us since it justifies the 
use of Bouleau-Hirsch's criterion (e.g. |24[ Section 2.1.2]) for establishing absolute 
continuity of F (cf. proof of theorem^). 

(2) Although not relevant to the sequel of this paper, it is interesting to compare 
Cy^-regularity with the Kusuoka- Stroock Sobolev spaces. Following |14[ I16j one de- 
fines D 1,p as the space of random-variables F which are (i) ray-absolutely-continuous 
(RAG) in the sense that for every h£TL there is an absolutely continuous version 
of the process {F (u + th) : t G M}; (ii) stochastically Gateaux differentiable (SGD) 
in the sense that there exists an TL-valued r.v. DF such that for every h G TL, 

(3.1) (F (uj +th) — F (w)) jt — » l\DF, as t -> 

probability with respect to p; and (Hi) such that F G L p and DF G L p (TL). 
From Sugita [27] it is known that D 1,p = D 1,p , at least for p G (l,oo). Since C^- 
regularity is a local property it has nothing to say about the integrability property 
(Hi) but it does imply a fortiori the regularity properties (i) and (ii). Indeed, (i) 
is trivially satisfied (without the need of h- dependent modifications!). As for (ii), 
DF is given by the Frechet differential DF (lu) of h G TL t— * F (u> + h) and the 
convergence H3.1\) holds not only in probability but p-almost surely. 
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Proposition 3. Let p < 3/2. For fixed t>0, the R e -valued random 



Me 



U: 



X(w) 



iVa) 



is continuously Ti-differentiable. 



Proof. Choose p > 2p such that l/p + 1/p > 1. We may assume that X (uj) has 
been defined so that X (uj) is a geometric p-rough path for every uj £ W. Let us also 
recall for h £ TL C C ,p_var , the translation T^X (uj) can be written (for uj fixed!) in 
terms of X (u>) and cross-integrals between -n\ (Xo,.) =: X £ C p ~ var and h. (These 
integrals are well-defined Young- integrals.) Thanks to the definition of X (ui) as 
the limit in probability of piecewise linear approximations to X and its iterated 
integrals (cf. [B]) and basic continuity properties of Young integrals we see that the 
event 



(3.2) 



{uj : X (uj + h) = T h X (uj) for all h £ TL} 



has probability one. We show that h £ TL ^ U. 



X(uj+h) 
t<-0 



(yo) is continuously Frechet 



differentiable for every uj in the above set of full measure. By basic facts of Frechet 
theory, we must show (a) Gateaux differentiability and (b) continuity of the Gateaux 
differential. 

Ad (a): Using X (uj + g + h) = T^X^X (uj) for g,h £ TL it suffices to show Gateaux 
differentiability of U i 



X(w+-) 



(yo) at £ TL. For fixed t, define 
Z i>s = J t l s (Vi (u£_o)) . 



Note that s 

i, 



Zi s is of finite p-variation. We have, with implicit summation over 



(yo) 



{Vi(uZ_ ))dh\ 



Z,dh' 



< 



< 



Z\ 



Z(Q)\) x \h\ 



\p—var 



Z \p-var + \ Z (°)l) >< \ h \n 



Hence, the linear map DU t 
component is an element of TL 

d 



o (y ) : h h-> D h U]t (y ) £ 
We just showed that 



is bounded and each 



de 



rr'T E ),X(u) 

u t^0 



(yo) 



DufSt ] (yo),h 



and hence 



a TT x(u+eh) , , 



H 



H 



emphasizing again that X (uj + h) = T^X (uj) almost surely for all h £ TL simultane- 
ously. Repeating the argument with T g X (uj) — ~K(uj + g) shows that the Gateaux 

differential of t/ t ^o +/) at g £ TL is given by 



DU, 



X(^+g) _ 



T s X(w) 



(b) It remains to be seen that g £ TL i— » DU, 



t^o 

T g X(w) 



£ L(TL,M. e ), the space of linear 



bounded maps equipped with operator norm, is continuous. To this end, assume 
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g in H (and hence in C p " var ). Continuity properties of the Young integral 
imply continuity of the translation operator viewed as map h 6 C p ~ val ^ T^X (cj) 
as p- rough path (see [5T]) and so 

in p-variation rough path metric. To point here is that 

x Jf_. and J£_. (Vj (U^o)) G C* p " var 

depends continuously on x with respect to p- variation rough path metric: using the 
fact that J£_. and U^_ both satisfy rough differential equations driven by x this 
is just a consequence of Lyons' limit theorem (the universal limit theorem of rough 
path theory). We apply this with x = X (u>) where u) remains a fixed element in 
(pT2"|) . It follows that 



s«-o 



DU, 



= sup 

op h:\h\ H =l 



DhU t +_ 



and defining Zf (s) = jj}^^ (Vt (u^o^jj, and similarly Zf n (s), the same 
reasoning as in part (a) leads to the estimate 



DU 



DU 



T 9 X(u 



< c 



op 



(V- 



Z 9 | 



I p-var 



+ (0) -Z 9 (0)|) . 



From the explanations just given this tends to zero as n — > oo which establishes 
continuity of the Gateaux differential, as required, and the proof is finished. □ 

Definition 2. gHHHIH] Given a continuously H- differential^ r.v. F = (F 1 , ...,F e 
W — > R e the Malliavin covariance matrix as the random matrix given by 

a(u) := ((DF\DF%). . =lt e eR**°. 

We call F weakly non- degenerate z/det (a) ^ almost surely. 

We now give an integral representation of the Malliavin covariance matrix of 
Y t = U^ ] (y ), the solution to the RDE driven by X (w), the lift of (X 1 , ...,X d ), 
along vector fields (Vi, ...,Vd), in terms of 2D Young integrals [3T | l29 l 18]. 

Proposition 4. Let a t = (J^DY t \ DY t ^ : i,j = 1, ej denote the Malliavin 

covariance matrix of Y t = uf^ (y ), the RDE solution ofdY = V (Y) dX (ui). In 
the notation of section [2.S\ we have 



DY t \ DYj 



d 



k=l 



JO 



J^ s (V k (Y s )) ® J t * s , (V k (Y s ,)) dR^ ( s , s ') 



Proof. Let [h { n ] : nj be an ONB of H {k) . It follows that (} 
is an ONB ofH = ® d l=1 U {k) where we identify 



(fc) 



1,2, ...;fc = 1,. 



V 



€ H 



J 



^^loc ls enou gh to guarantee the existence of an W-valued derivative. 



NON-DEGENERACY OF WIENER. FUNCTIONALS ARISING FROM RDES 



and similarly for k = 2, d. From Parseval's identity, 
a t = ((DY*,DY t \ . 



^{DYuhM) H ®(DY t ,hg> i 

n,k 

E E f ^ ^)) ® f ^ ^ dh n-s 

k n J JO 

E f f J^ s {v k {Y s ))® j* s , (v fc (y y ))dflW (*,*'). 

,, Jo Jo 



k 

For the last step we used that 

f(s)g(t) dR(s,t) 



E / f dhn / 9dhn = I \ 
„ Jo Jo Jo Jo 



whenever f — f (t, u>) and g are such that the integrals are a.s. well-defined Young- 
integrals. The proof is a consequence of R (s, t) — E (X s X t ) and the L 2 -expansion 
of the Gaussian process X, 

X(t) = Y,^(hn)h n (t) 

n 

where £ (h n ) form an IID family of standard Gaussians. □ 



Two special cases are worth considering: in the case of Brownian motion dR (s, s') 
is a Dirac measure on the diagonal {s — s'}) and the double integral reduces to a 
(well-known) single integral expression; in the case of fractional Brownian mo- 
tion with H > 1/2 it suffices to take the d 2 / (dsdt) derivative of Rn(s,t) — 

f f 2H + S 2H _ | t _ s |2H^j j 2 tQ gee that 

dR H (s,s') ~ \t - s\ 2H ~ 2 dsdt 

which is integrable iff 2H — 2 > — 1 or H > 1/2. (The resulting double- integral rep- 
resentation of the Malliavin covariance is also well-known and appears, for instance, 
in [221C11IT].) 

4. Existence of a density for Gaussian RDEs 

We remain in the framework of the previous sections, Y t (to) = U^2q (yo) de- 
notes the (random) RDE solution driven a Gaussian rough path X, the natu- 
ral lift of a continuous, centered Gaussian process with independent components 
X = (X\...,X d ) started at zero. Under the standing assumption of finite 
/9-variation of the covariance, p < 3/2, we know that X (cj) is a.s. a geometric 
p- rough path for p £ (2p, 3). Recall that this means that X can be viewed as path 
in G 2 (R d ), the step-2 nilpotent group over R d , of finite p-variation relative to the 
Carnot-Caratheodory metric on G 2 (R d ) . 

Condition 1. Ellipticity assumption on the vector fields: The vector fields 
V\, Vd span the tangent space at yo. 
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Condition 2. Non- degeneracy of the Gaussian process on [0, T\: FixT > 0. 
We assume that for any smooth f = (/i, fd) '■ [0, T] — > K.'' 




Note that non-degeneracy on [0, T] implies non-degeneracy on [0, t] for any t £ 
(0,T]. It is instructive to see how this condition rules out the Brownian bridge 
returning to the origin at time T or earlier; a Brownian bridge which returns to 
zero after time T is allowed. The following lemma contains a few ramifications 
concerning condition [2] Since Ti = (B d k=1 H^ there is no loss in generality in 
assuming d = 1. 

Lemma 1. (i) The requirement that f is smooth above can be relaxed to f £ C p ~ var 
for p > 2p small enough. 

(ii) The requirement that J fdh = OVTi £ Ti can be relaxed to the the quantifier 'for 
all h in some orthonormal basis of Ti" . 

(Hi) The non- degeneracy condition^ is equivalent to saying that for all smooth 
f 7^ 0, the zero-mean Gaussian random variable J Q fdX (which exists as Young 
integral or via integration-by-parts) has positive definite variance.) 
(iv) The non-degeneracy condition^ is equivalent to saying that for all times < 
ii < ... < t n < T the covariance matrix of (X tl , X tn ), that is, 

is (strictly) positive definite. 



Proof, (i) Continuity properties of the Young integral, noting that 1/p + l/p > 1 
for p > 2p small enough, and weak compactness of the unit ball in Ti.. (ii) Any 
h e H can be written as limit in Tt of h^ = 53fc=i {hk, h) hk as n — > oo when (hk) 
is an ONB for TL. From Proposition [2] it follows that h is also the limit of foW in 
p-variation topology. Conclude with continuity of the Young integral, (iii) Every 
element h £7i can be written as h t = E (ZX t ) with Z £ £ (H). By taking L 2 -limits 
one easily justifies the formal computation 



dht = 



(x t z) 



dt 



fdh = E 



fdX Z 



and our condition is equivalent to saying that 



f fdX±£(H) 
Jo 



On the other hand, it is clear that J fdX itself is an element of the 1st Wiener Ito 
Chaos £ (Ti.) and so must be in I? . In other words, saying that J fdh — for all 
h £ Ti says precisely that 

2 











Var 


f fdX 




f fdX 




Jo 




Jo 



= 0. 



L- 



Conversely, assume that 



Var 



fdX 
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Then J fdX = with probability 1 and by the Cameron-Martin theorem, for any 

hen, 

fd(X + h) = Q 

with probability one which implies J fdh = for all h E H. (iv) Suffices to note 
that 

[ T [ T f (*) f (t) dR (s, t) = V f f (s) dh, 
Jo Jo „ Jo 



For the last step we used that 



y~] / fdh n / gdh n = I / 
„ Jo Jo Jo Jo 



f(s)g(t)dR(s,t) 



which is a consequence of R(s,t) — E,(X s X t ) and the L 2 -expansion X (t) — 
J2 n £ C 1 ") hn (t). (v) Left to the reader, (vi) Obvious. □ 

Remark 2. The variance of J fdX can written as 2D Young integral, 

f s f t dR(s,t). 

[0,T] 2 

To put the following result in context, recall that the covariance of Brownian 
motion, (s,t) i— > min(s,£), has finite 1-variation with p = 1. For fractional Brown- 
ian motion with Hurst parameter H one can take p = 1/2H. The following result 
then applies to fractional Brownian driving signals with H > 1/3. 

Theorem 1. Let X be natural lift of a continuous, centered Gaussian process with 
independent components X — (X 1 , X d ^j , with finite p-variation of the covari- 
ance, p < 3/2 and non-degenerate in the sense of condition^ Let V — (V±, Vd) 
be a collection of Lip 3 -vector fields on R e which satisfy the ellipticity condition^ 
Then the solution to the (random) RDE 

dY = V{Y)dX, Y (0) = g G M e 

admits a density at all times t € (0, T] with respect to Lebesgue measure on R e . 

Proof. Fix t G (0, T]. From Proposition [3] we know that U^^yo — Y t is continu- 
ously 7i-differentiable. By a well-known criterion due to Bouleau- Hirsdfl the proof 
is reduced to show a.s. invertibility of the Malliavin covariance matrix 

<j t = ((dY^DyA ) £ M exe . 

Assume there exists a (random) vector u£l e which annihilates the quadratic form 
a t . TherQ 

2 







T 

V 0~tV 



J2 ViDY t 



and so v T DY t = ^ v^DY} eOeH. 



By propositions [T] and [3l 

(4.1) VheH: v T D h Y t = f V v T j*_ a {V 3 (Y s )) dh? s = 

Jo , =1 



ft d 



^Combine the result of [Nualart, 4.1.3] and [Nualart, section 2]. 
^Upper T denotes the transpose of a vector or matrix. 
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where the last integral makes sense as Young integral since the (continuous) inte- 
grand has finite p-variation regularity. Noting that the non-degeneracy condition 
on [0, T] implies the same non-degeneracy condition on [0, t] we see that the inte- 
grand in (|4.1[) must be zero on [0, t] and evaluation at time shows that for all 
j = 1, ...,d, 

v T J?io (V, M) = 0. 
It follows that the vector v T J£_ is orthogonal to Vj (yo) , j = 1, ...,d and hence 
zero. Since J^_ is invertible we see that v = 0. The proof is finished. □ 

The reader may be curious to hear about smoothness in this context. Adapting 
standard arguments would require LP (f2) estimates on the Jacobian of the flow 
J^o' '. Using the fact that it satisfies a linear RDE, dj£_ to (y Q ) = c?M x • j£_ t (yo), 
with cZM x = V (Y) dX one can see that 

(4-2) log 1^(2/0)1=0(11X11^.). 

(This estimate appears in (55) for p < 2 but can be seen [ISHHj to hold for all p > 1. 
We believe it to be optimal.) Using the Gauss tail of the homogenous p-variation 
norm of Gaussian rough paths (see [7J [5] ) we see that L 9 -estimates for all q < 00 
hold true when p < 2 and this underlies to density results of [26l Q] . On the other 
hand, for p > 2 one cannot obtain L 9 -estimates from (|4.2|) and further probabilistic 
input will be needed. 
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